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ROLE  OF  CARBON  DIOXIDE  IN  COOLING  PLANETARY  THERMOSPHERES 

Ramesli  D.  S harm  a 
Geophysics  Laborator}' 

Peter  P.  Wintersteiner 
ARCON  Corp. 


Abstract.  A  new  value  of  the  rate  coefficient  for  the  de¬ 
activation  of  the  bending  mode  of  carbon  dioxide  by  ato¬ 
mic  oxygen  at  low  temperatures  is  derived  from  the  obser¬ 
vation  of  15  /im  emission  from  the  atmosphere  of  the 
Earth.  This  new  value  gives  a  cooling  rate  for  the  lower 
thermosphere  that  is  two  to  three  times  the  rate  previous¬ 
ly  calculated,  and  it  may  resolve  a  long-standing  problem 
in  the  Mars-Venus  aeronomy. 

Introduction 

The  role  of  carbon  dioxide  in  cooling  the  thermospheres 
of  Earth,  Venus,  and  Mars  is  well  recognized  [Gordiets  et 
al,  19S2,  Dickinson,  19S4,  Dickinson  et  al,  1987;  Dickinson 
and  Bougher,  1986,  Bougher  and  Dickinson,  198S].  During 
a  collision  of  C02  with  atomic  and  molecular  species, 
some  of  the  translational  energy  of  relative  motion  (heat) 
is  converted  into  vibrational  energy  of  the  lowest  lying 
mode,  the  bending  mode,  i/2.  This  vibrational  energy  is 
subsequently  radiated  at  15  /<m.  Many  of  these  photons 
escape  to  space,  cooling  the  atmosphere.  A  critical 
parameter  in  this  process  is  the  rate  coefficient  foi  the  de¬ 
activation  of  the  bending  mode  by  atomic  oxygen.  Tayloi 
[1974]  gives  its  value  as  2.32x  10-9exp(  76.75/T1/3) 
cm3/mol-s.  This  expression  has  the  standard  Landau- 
Teller  form  and  is  obtained  by  extrapolation  from  high 
temperature  shock  tube  experiments.  There  arc  no  low- 
temperature  measurements,  however,  and  the  late 
constant  is  considered  to  be  very  uncertain. 

We  suggest  that  a  revised  value, 

k0=  3.5(±1.8)xlO_13'[T  + 

2.32  x  10"9exp(-76.75/T,/3)  (1) 

better  leprcsents  the  overall  temperature  behavior  and 
should  be  used  for  atmospheric  studies.  This  revised  value 
is  larger  than  that  given  by  Taylor  by  a  factor  of  1070  at 
200  °K  andbv  a  factor  of  54  at  500  °K.  At  300  °K  it 
equals  6x10  13cm3/mol-s.  Sharma  and  Nadile  [1981]  in  a 
preliminary  study  proposed  5xl0~13  for  this  rate  constant 
at  room  temperature.  Values  close  to  this  have  since  been 
used  in  the  literature  [Gordiets  et  al,  1982;  Roble  et  al, 
19S7;  Dickinson  and  Bougher,  1986;  Kumer  and  James, 
1982],  The  value  being  proposed  in  this  article  is  about  12 
times  larger  than  the  preliminary  value  and  is  closer  to 
results  from  the  expression  2xl0"13>jT  given  by  Sharma 
[198S]  for  low  temperatures. 

Enhanced  emission  resulting  from  a  large  upward  revi¬ 
sion  of  this  rate  coefficient  affects  the  modeling  of  the 
mesosphere  and  lower  thermosphere  of  the  Earth.  The 
impact  on  the  understanding  of  the  energy  budget  of 
Venus  may,  however,  be  dramatic.  The  upper  atmos¬ 
pheres  of  Venus  and  Mars  are  composed  mainly  of  C02, 
with  small  admixtures  of  0,  CO,  and  02/N2  in  decreasing 
proportion  [Chamberlain  and  Hunten,  1987].  While  the 
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thermospheric  temperature  structure  is  reasonably  well 
understood  for  Mars,  the  corresponding  understanding  for 
Venus,  which  has  a  similar  atmosphere,  is  lacking  [Fox 
and  Bougher,  in  press,  1990].  (The  major  compositional 
difference  is  the  mixing  ratio  of  0,  ~  1%  at  the  Martian 
homopause  [Neir  and  McElroy,  1977]  compared  to  about 
6-S%  on  Venus  [Hedin  et  al,  1983].)  Fox  [1938]  has 
calculated  large  EUV  heating  efficiencies  (16-25%)  for 
Venus  which,  in  turn,  result  in  predictions  of  the 
exospheric  temperature  that  are  much  too  high. 

Recently  Keating  and  Bougher  [in  press,  1990]  have 
pointed  out  that  the  response  of  the  Venus  thermosphere 
to  short-term  solar  variations  can  only  be  explained  by 
very  strong  cooling  in  the  15  /im  bands,  and  not  by  very 
large  eddy  cooling  [cf.  Hollenbach  et  al,  1985].  They  also 
point  out  that  the  heating  efficiency  calculated  by  Fox 
[19SS]  together  with  the  late  coefficient  kG  proposed  here 
yield  a  leasonable  model  for  the  Venetian  thermosphere. 

SPectral  Infrared  Rocket  Experiment 

SPIRE  was  launched  from  the  University  of  Alaska’s 
Poker  Flat  Research  Range  shortly  before  dawn  on 
September  2S,  1977,  into  a  quiescent  non-auroral  atmos¬ 
phere  [Stair  et  al,  19S5).  Foi  about  five  minutes  it  mea¬ 
sured  limb  emission  in  the  spectial  range  1.4-16.5  tan.  It 
made  twelve  spatial  scans.  Scans  1-7  were  termmatoi 
scans  for  which  the  lines  of  sight  (LOS)  were  partially 
sunlit,  scan  8  was  a  night-time  scan,  and  scans  9-12  had 
LOS  completely  sunlit. 

A  plot  of  radiance  in  the  r2  band  of  C02  as  a  function 
of  tangent  height  shows  several  features,  of  which  two  are 
worthy  of  special  notice.  (1)  Within  the  scatter  of  the 
data,  the  day  and  night  measurements  do  not  show  any 
systematic  differences.  In  other  words  the  15  /an  emis¬ 
sion,  early  in  the  morning,  is  independent  of  whether  the 
LOS  arc  dark  or  sunlit.  (2)  The  band  radiance  remains 
nearly  constant  as  the  tangent  height  increases  from 
about  94  to  110  km.  This  is  important,  because  according 
to  all  models  of  the  atmosphere  the  C02  density  should 
decrease  by  a  factor  of  at  least  thirty  in  this  altitude 
range.  Recalling  the  fact  that  the  radiance  observed  on 
any  tangent  path  originates  primarily  at  the  tangent 
altitude,  and  ignoring  small  effects  due  to  self- absorption, 
we  conclude  that  the  excited  state  density  stays  almost 
constant,  or  the  ratio  of  the  densities  of  the  upper  radia¬ 
ting  state  (01  ’0)  to  the  lower  radiating  state  (00°0)  of 
C02  increases  by  a  factor  of  thirty  or  more,  as  we  go  from 
94  to  110  km. 

Model  of  the  15  /an  Emission  from  C02 

We  recall  the  observation  that  solar  radiation,  early  in 
the  morning,  does  not  play  an  important  role  in  the 
emission  process  at  15  /an.  Consider  the  ground  and  first 
excited  states,  00°0  and  01*0.  The  molecules  in  the  upper 
state  are  produced  from  those  in  the  lower  state  by  colli¬ 
sions  with  N2,  02,  and  0  and  by  absorption  of  the  upwel- 
ling  radiation  from  the  atmosphere  (earthshine).  The 
upper-state  molecules  N„  are  converted  to  the  lower-state 
molecules  N(  by  collisions  with  N2,  02,  and  O  and  by 
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spontaneous  emission.  The  steady-state  equation  for  the 
two-level  system  may  then  be  written  as 

Nu/N,  =  (J  +  kK)/(A  +  k)  (2) 

where  J  is  the  pumping  rate  due  to  earthshine,  A  is  the 
Einstein  coefficient  for  spontaneous  emission,  k  = 
£kM[M)  and  kK  are  the  collisional  de-excitation  and 
excitation  rates  of  the  upper  state,  K  =  2  exp(-960/T)  is 
the  equilibrium  constant  (960  °K  being  the  energy 
separation  of  the  two  levels),  T  is  the  translational 
temperature,  and  the  factor  of  two  in  the  expression  for  K 
arises  from  the  ratio  of  statistical  weights  of  the  upper 
and  lower  states. 

The  rate  of  collisional  deactivation  of  Ol'O  by  N2  and 
02  has  been  measured  in  the  temperature  range  of 
interest  [Allen  et  al,  19S0;  Lunt  et  al,  1985;  Simpson  et  al, 
1977;  Taine  et  al,  197S;  Taine  and  Lepoutre,  1979).  The 
contribution  of  atomic  oxygen  has  already  been  men¬ 
tioned.  It  is  the  computation  of  J  that  forms  the  heart  of 
the  calculation.  As  the  upwelling  radiation  leaves  the 
lower  atmosphere,  the  photons  at  the  line  centers  of 
strong  bands,  like  the  one  under  discussion,  are  rapidly 
absorbed  or  bleached  out.  A  large  fraction  of  the  radiation 
contributing  to  J  comes  from  the  line  wings.  The  radia¬ 
tive  transfer  calculation  must  be  carried  out  carefully, 
taking  into  account  changes  in  the  lineshape  and  the 
absorption  coefficient  as  the  temperature  and  density 
change  with  altitude.  We  wrote  an  accurate  line-by-line 
code,  called  RAD,  for  the  purpose  of  calculating  J 
(Wintersteiner  et  al,  19SS).  The  code  assumes  a  one¬ 
dimensional  plane-parallel  atmosphere  and  that  the 
spectral  lines  do  not  overlap.  It  equates  the  rotational 
temperature  with  the  translational  temperature  and 
assumes  that  an  absorbed  photon  may  be  emitted  any¬ 
where  in  the  line  with  a  probability  determined  by  the 
shape  of  the  line.  RAD  calculates  the  vibrational  tempera¬ 
ture  from  the  logarithm  of  the  right  hand  side  of  Eq.  (2), 
iteratively. 

Model  Atmosphere 

Since  the  rocket  was  at  approximately  2S5  km  while 
the  measurements  were  being  made,  in-situ  measurements 
of  the  atmospheric  parameters  could  not  be  carried  out. 
We  used  the  MSIS-S6  thermospheric  model  [Hedin,  19871 
to  give  the  temperature,  densities  of  N2,  02.  and  0,  and 
the  total  density  under  the  SPIRE  conditions.  Below  the 
85-km  lower  boundary  of  MSIS,  the  sub-arctic  winter 
atmosphere  given  in  FASCODE  [Anderson  et  al,  19S6] 
was  used,  with  suitable  interpolation  near  the  boundary. 
We  were  able  to  correctly  reproduce  the  rotational  profile 
of  the  15  pm  band  using  the  temperature  thus  obtained, 
giving  some  credibility  to  this  procedure.  For  [0]  we  used 
a  simple  extrapolation  of  the  MSIS  profile  below  85  km. 
Figure  1  shows  (0)  as  a  function  of  altitude. 

To  estimate  [C02]  we  determined  mixing  ratios  from 
all  thermospheric  measurements  reported  in  the  litera¬ 
ture.  Above  110  km,  the  measurements  are  in  general 
agreement,  but  in  the  altitude  range  80-110  km  there  is  a 
disparity  between  the  results  of  rocket  experiments 
[Trinks  and  Fricke,  1978;  Trinks  et  al,  1978;  Offermann  et 
al,  1981;  Offermann  and  Grossmann,  1973]  and  satellite 
experiments  [Vercheval  et  al,  1986;  Girard  et  al,  19S8; 
Beer  et  al,  1987),  with  the  former  measurements  consis¬ 
tently  showing  up  to  50%  more  C02  than  the  latter.  We 
represented  the  data  with  two  smoothed  mixing  ratio 
profiles  and  ran  the  model  using  C02  densities,  shown  in 
Figure  1,  that  were  derived  from  each  of  them.  Curve  D 
refers  to  the  profile  derived  from  the  rocket  data,  and 
curve  E  to  that  from  satellite  data. 
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Fig.  1.  Atomic  oxygen  and  carbon  dioxide  densities  as  a 
function  of  altitude.  The  [0]  profile  is  from  MSIS-86; 
[C02]  curves  D  (solid  line)  and  E  (dashed  line)  are  approx¬ 
imations  to  rocket  and  satellite  measurements,  respec¬ 
tively. 


Comparison  of  Measured  and  Calculated  Radiances 

After  computing  the  vibrational  temperatures,  we  used 
a  line-by-line  radiative  transport  code,  NLTE  [Wmter- 
steiner  and  Sharma,  1985],  to  compute  the  limb  radiance 
in  each  rotational  line  of  the  Ol'O— 00°0  vibrational  bands 
of  the  four  most  important  isotopes  of  C02.  RAD  and 
NLTE  were  also  used  to  compute  the  limb  radiance  from 
four  hot  bands  (0330-02*0,  10°0—0110,  02°0-0l'0, 
0220— Ol'O)  of  the  most  abundant  isotope,  l2C1602  and 
fiom  the  0220— Ol'O  band  of  13C,602.  At  85  km  tangent 
height  the  major  and  minor  isotopes  make  equal  contiibu- 
tions  to  the  radiance  in  the  13-16.5  ,im  wavelength  inter¬ 
val  because  of  greater  self-absorption  of  radiation  from 
the  most  abundant  isotope.  Above  105  km,  however,  only 
the  Ol'O— 00°0  transition  of  the  most  abundant,  isotope 
makes  a  significant  contiibution  to  the  radiance  at  15  /«m. 

We  ran  the  model  using  the  two  C02  profiles  and 
compared  the  predicted  limb  ladiance  with  the  SPIRE 
dataset.  The  94-110  km  featuie  discussed  earlier  is  due  to 
the  efficiency  with  which  O-atom  collisions  excite  the 
bending  mode,  this  is  exactly  the  region  where  [O]  and  k0 
increase  rapidly  with  altitude.  It  is  impossible  to  repro¬ 
duce  this  feature  with  accepted  values  of  the  rate  con¬ 
stant,  so  we  let  k0  (or,  more  precisely,  the  coefficient  of 
\Tr  in  Eq.  (1))  be  a  free  parameter.  Figures  2  and  3  com¬ 
pare  the  measured  radiance  with  three  radiance  profiles, 
calculated  with  different  0-atom  rate  constants  having 
the  form  A  x  10“l3'jT  +  2.32  x  10~9exp(-76.75/T1/3).  The 
best  overall  description  of  the  observations  is  obtained 
with  a  value  A  -3.5. 

Sensitivity  to  Model  Input 

In  our  model,  [0]  and  k0  appear  only  as  a  product,  so 
any  error  in  the  former  is  immediately  reflected  in  the 
result  we  obtain  for  the  latter.  Although  [0]  is  quite 
variable,  the  standard  deviations  between  the  MSIS-86 
model  and  lower-thermosphere  data  reportedly  range  up 
to  about  30%  [Hedin,  19SS).  If  the  MSIS  predictions  were 
known  to  b"  in  error  by  substantially  more  than  this,  our 
error  estimate  for  ko  would  have  to  be  correspondingly 
revised. 
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SPIRE  BAND  RADIANCE,  13-16.5  fJtm 


Fig.  2.  Limb  radiance  in  the  13.0-16.5  /im  band  as  a  func¬ 
tion  of  tangent  height.  Asterisks  denote  the  SPIRE  mea¬ 
surements.  The  calculated  radiances  use  C02  profile  D. 
The  parameterization  for  kG  is  Ax10~,:MT  +  2.32  xlO-9 
exp(--76.75/T'^3),  and  the  values  used  for  A  are  indicated. 


SPIRE  BAND  RADIANCE,  13-16.5  /zm 


Fig.  3  As  Figure  2,  but  the  radiances  were  calculated 
using  C02  profile  E. 


Above  110  km,  the  limb  radiance  is  proportional  to  the 
C02  column  density  on  the  viewing  path.  In  this  altitude 
range  the  mixing  ratio  dataset  has  a  standard  deviation  of 
about  25%  about  the  mean.  Test  runs  show  that  a  25% 
change  in  [C02]  is  compensated  by  a  change  of  approxi¬ 
mately  36%  in  kG.  For  lower  tangent  paths,  even  though 
the  uncertainty  in  [C02]  is  greater,  self-absorption  sharply 
diminishes  the  sensitivity  of  the  limb  radiance  to  the 
column  density. 

The  limb  radiance  in  the  thermosphere  is  also  sensitive 
to  the  temperature  profile,  particularly  to  the  mesopausc 
temperature.  We  ran  our  model  with  two  FASCODE 
atmospheres  (Anderson  et  al,  19S6]  having  lower  and 
higher  mesopause  temperatures.  Radiance  results  ob¬ 
tained  using  the  midsummer  profile,  with  a  cold  meso¬ 
pause,  are  consistent  with  the  experimental  data  when  kG 
is  given  by  Eq.  (1)  However,  results  obtained  with  the 
midwinter  profile,  with  a  warmer  mesopause,  do  not 
match  the  experiment  and  cannot  be  reconciled  with  it 


by  adjusting  kG  downward.  This  indicates  that  the 
temperature  profile  is  at  least  approximately  correct,  and 
that  kc  cannot  be  greatly  in  error  on  this  account. 

The  overall  effect  of  the  major  uncertainties  in  model 
input  leads  to  an  error  of  xiproxim.'CCiy  50%  in  the  rate 
coefficient,  as  given  in  Eq.  (1). 

Cooling  Rates 

Figure  4  gives  the  cooling  rate  due  to  15  pm  emission 
as  a  function  of  altitude  for  both  C02  profiles.  Below  120 
km  this  is  the  principal  cooling  mechanism.  Since  it  is 
very  sensitive  to  k^,  these  results  represent  a  substantial 
increase  over  cooling  rates  previously  reported  [cf. 
Dickinson,  19S4J.  This  has  important  implications  for  the 
transport  of  heat  into  and  within  the  lower  thermosphere, 
and  for  its  thermal  balance  as  well. 


15  fj.m  COOLING  RATES 


Fig.  4.  Cooling  rates  due  to  15  pm  emission  from  the 
Earth’s  atmosphere,  as  a  function  of  altitude,  for  both 
C02  profiles. 


Summary 

Based  upon  observations  of  15  /an  emission  from  the 
atmosphere  of  the  Earth,  we  have  deduced  a  new  value 
for  the  rate  of  deactivation  of  the  bending  mode  of  carbon 
dioxide  by  atomic  oxygen.  Using  this  rate  coefficient  we 
calculate  cooling  rates  for  the  terrestrial  thermosphere 
that  arc  larger  than  generally  recognized.  This  will  also 
give  larger  cooling  rates  for  the  Venus  thermosphere,  and 
possibly  resolve  the  long-standing  problem  regarding  its 
temperature  structure. 

Because  of  its  extremely  important  role  in  cooling  pla¬ 
netary  thermospheres,  it  is  of  great  interest  to  obtain 
independent  confirmation  of  these  low- temperature 
results  in  the  form  of  laboratory  measurements. 

Acknowledgement.  The  authors  are  grateful  to  Prof. 
Jane  L.  Fox  for  discussions  about  the  Mars  and  Venus 
aeronomy. 
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